Surface instability and isotopic impurities in quantum solids 
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In this paper we employ a self-consistent harmonic approximation to investigate surface melting 
and local melting close to quantum impurities in quantum solids. We show that surface melting 
can occur at temperatures much lower than the critical temperature Tc of the solid phase instability 
in the bulk. Similar effects are driven by the presence of an isotope substitution. In this latter 
case, we show that stronger local lattice fluctuations, induced by a lighter isotope atom, can induce 
local melting of the host bulk phase. Experimental consequences and the possible relevance in solid 
helium are discussed. 



I. INTRODUCTION 

Although melting is a very common phenomenon in 
nature, the debate about its microscopic mechanism is 
still openii>^i^ The first empirical theory was advanced 
by Lindemann.— According to this view, melting occurs 
when the ratio between the root mean square (rms) 
Mrms = \/ (u^) of the thermally activated lattice fluctu- 
ations and the lattice constant a exceeds a phenomeno- 
logical threshold Mrms/a <^ 0.16 which is roughly material 
independent^'^. In spite of its several flaws (melting is de- 
scribed in terms of the properties of only the solid phase; 
no cooperative process and no role of defects are con- 
sidered. . . ), this simple criterion seems to work reason- 
ably well for a variety of materials.^ Lindemann criterion 
has been recently found to apply as well at a local level 
around crystal defects.— 1^ This large range of validity of 
the Lindemann criterion suggests thus that a microscopic 
mechanism is actually operative. 

The most simple (and employed) model to account 
for the Lindemann phenomenology is the self-consistent 
harmonic approximation (SCHA). This maps an anhar- 
monic phonon model in a harmonic one. Anharmonic- 
ity is, then, taken into account, at a mean-field level, 
through a Debye- Waller-like term which is evaluated self- 
consistently. The breakdown of this approach is inter- 
preted as a signal of solid phase instability, and hence 
related to melting. One of the strength of this theory is 
that it predicts, in contrast with the Born criterion but 
in agreement with the experimental observation, a par- 
tial but not total softening of the elastic constants of the 
bulk. 

The SCHA represents moreover an efficient tool to un- 
derstand in a qualitative way the phenomenon of the sur- 
face melting (SM), as first proposed by Pietronero and 
Tosatti (PT).^" In this context, the physical mechanism 
underlying the surface melting is quite simple: atoms 
close to the surface have larger lattice fluctuations due 
to the reduced number of nearest neighbor sites, and 
the SCHA breaks down consequently at smaller temper- 
atures than in the bulk. It is clear that this simple theory 
does not represent an exhaustive description of the sur- 



face melting phenomenology, which should include rough- 
ening, preroughcning, partial wetting, the role of "crys- 
tallinity" etc.^^-^^ In addition it should be stressed that 
the SCHA does not determine directly the melting point 
but rather the instability of the solid phase which is pre- 
vented by the melting process itself. ^'^ In this perspective 
this criterion should not be employed at a quantitative 
level. Nevertheless, since the solid phase instability and 
the actual melting process are usually related to each 
other, the PT theory provides a simple and useful way 
to get information about the tendency of a system to- 
wards melting and surface melting and their dependence 
on microscopic parameters. 

In this paper we generalize the results of the PT ap- 
proach in the case of quantum solids. The Lindemann 
criterion in the quantum solid is shown to be twice as 
large as the one in the classical limit, in agreement with 
experimental reports We show a phase diagram for 
both the bulk and surface melting cases and we investi- 
gate also the local melting due to an isotopic substitution. 
The temperature dependence of the lattice fluctuations 
for the different classical/quantum regimes is evaluated 
and also the profile of the lattice fluctuations as function 
of the distance from the surface or the isotopic impurities. 
The paper is organized as follows: in Sect. |TT1 we review 
the approach of Pietronero- Tosatti for classical solids; in 
mil we generalize the PT approach to investigate bulk 
properties in quantum solids; surface melting and solid 
phase instability close to a quantum isotope impurity are 
analyzed respectively in Sect. IIVI and Sect. |Vl Finally, 
in Sect. IVIi we discuss our results and draw some flnal 
conclusions. 



II. SCHA AND SOLID PHASE INSTABILITY IN 
BULK AND ON SURFACES 

Let us consider for simplicity a one-dimensional chain 
of atoms. At the harmonic level, we can write the equa- 
tions of motion for the lattice displacement u„: 
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where M is the atomic mass, n denotes the site in- 
dex. The constant forces kn^„i, at the harmonic level, 
are related to the inter-atom potential Vn.n' through the 
relation kn^n' — (?^Ki,ri'/^^nl{M„}=o- Writing the po- 
tential Vn,n' in terms of a Fourier expansion, Vn,n' — 



')], we have thus, at the harmonic 



level, kn,n' =ko = -J2q 

In the spirit of SCHA, anharmonic terms can be taken 
into account, by replacing the constant forces fc„_„', eval- 
uated at the lattice equilibrium, with their expectation 
value fc„,„' averaged over the lattice fluctuations. We 
have thus explicitly: 
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where in the last line we have neglected the cross terms 
and we have replaced the dependence on the momenta in 
the exponential with an effective parameter A. 

By inserting ([2|) in Eq. ([T]) and considering the motion 
of each atom as an Einstein oscillator we have: 
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where anharmonic effects are taken into account in the 
self-consistent renormalization of the elastic constants 
kn.n> ■ Note that fc„^„' depends on the expectation value 
of the quadratic lattice fluctuations on both sites n, n' . 
It follows that the atomic motion described in Eq. ([3]) 
is ruled by the lattice fluctuations of the lattice environ- 
ment. In a bulk system (u^) = {u^^i) = {'^^)) then 



force constant fco, in agreement with the observation of 
a material independent Lindemann criterion. 

Eq. ^ represents also the starting point to apply the 
SCHA to surface melting. In this case, one defines a local 
average lattice fluctuation (it^) which depends on the site 
index n. In the same spirit one can define a local elastic 
constant: 



^n,n — l,n+l ^" ^n,n — 1 
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We can write thus a set of recursive equations where the 
lattice fluctuations of the atom n depend on the lattice 
fluctuations of the n — 1 and 71 + 1 atoms. The recur- 
sion is truncated at the atom n — \ which represents the 
outer atom close to the free surface. This atom probes 
an effective harmonic potential smaller than the bulk, 
which increases its tendency towards melting. A numer- 
ical solution shows that the solid phase for the surface 
atoms becomes unstable at r™ = 0.271, 26 % smaller 
than the bulk value. The same theory permits to evalu- 
ate the profile of the lattice fiuctuations as function of 
the distance from the surface. These theoretical pre- 
dictions agree quite well with the profile of the lattice 
fluctuations close to defects (grain boundaries, disloca- 
tions, vacancies) in colloidal solidSfS Note that, although 
the temperature of surface melting is smaller than in the 
bulk, local lattice fluctuations of the outer atoms can be 
larger than the ones in the bulk, violating locally the 
Lindemann criterion. This is also in agreement with Rcf. 
[9^1 . For instance, for the outer atoms n = 1 one flnds 
ySM ^ i j^^ rpj^-g larger than the value in the 

bulk. 



kn,n' = k = ko cxp [-A(u^)] . 
and we get an unique self-consistent equation 
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where fee is the Boltzmann constant. In similar way. 



the SCHA phonon frequency is given hy ujq = y k / M = 

Wo exp[— A(m^)/2], where luq = \/k()/M is the bare 
phonon frequency at the purely harmonic level. It is 
convenient to rewrite Eq. ^ by introducing the dimen- 
sionless quantities y = A(w^), Td = \k^T/ko: 
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Eq. © has no solution for Td > t^^^^ ~ 1/c = 
0.368, which determines a critical temperature k^Tc — 
0.368fco/'^- At this value y{T^^^) = 1 and the maximum 
magnitude of the allowed lattice fluctuations above which 
the solid phase is unstable is (u^)™'^'' = 1/A. Note that 
^y2^max (iepgj^fjs neither on the atomic mass nor on the 



III. BULK PROPERTIES OF QUANTUM 
SOLIDS 



We generalize now the above theory to the case of 
quantum solids. In the following we shall assume a one- 
particle picture to be still valid, because of the smallness 
of the exchange terms in the solid phase (J™'^'^ '--^ 0.1 K 
in ^He, J™'*'^ ~ in "^He) with respect to the melt- 
ing temperatures Tm > 2 KJ^^ On the other hand a 
major role in our approach will be played by the quan- 
tum fluctuations which dominate at low temperature in 
the quantum regime. According to this perspective, the 
atomic motion of the atom n is described in terms of the 
SCHA Hamiltonian of the quantum oscillator: 
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where the self-consistent expression for the local potential 
fcn,n-i,n+i is reported in Eq. ([7]). 

We consider first the melting properties of bulk sys- 
tems {kn^n-i,n+i = 2fc). In this SCHA quantum model 
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the total amount of lattice fluctuations is now easily com- 
puted as: 



2Mujo 



(9) 



where n{x) = l/[e^ — 1] is the Bose factor and where we 

remind (Dq — \Jk/M and k is given by EqU In the classic 

limit k^T ^ Tiujo^ n{x) ~ 1/a; 3> 1, and we recover the 
classical result of Eq. ([5]). On the other hand, in the 
zero temperature limit, lattice fluctuations are due only 
to zero point quantum motion. In this case n{x) — and 
Eq.® reads: 

:exp[A(7/2)/2] , (10) 



which, introducing the variable tq — Xh/2^/k^M, can be 
written in the dimensionless form: 



y(TQ) = TQe2'(^«)/^ 



(11) 



Eq. (jlip represents the quantum generalization of Eq. 
([6]) where the instability of the solid phase is now trig- 
gered by the magnitude of the quantum lattice fluctua- 
tions. This occurs for tq > t^^"" = 2/e = 0.736. It is 
interesting to note that the breakdown of the solid phase 
driven by quantum fluctuations is not merely equivalent 
to the one related to the thermal motion. Indeed, for a 
quantum solid, we would predict a maximum magnitude 
of lattice fluctuations 2/(tq '^'^) = 2, two times larger than 
for classical solids. This behavior is indeed in agreement 
with the report of the Lindemann ratio Umis/a — 0.28 in 
helium solidsi^iiiii^ to compare with itmis/a — 0.16 for 
classical solids. 

We also consider now the general case where both ther- 
mal and quantum fluctuations are important. From Eq. 
([5]), after few straightforward passages, we get 



2/(t-q,Tci) = Tqe 
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Eq. (fT^ generalizes the stability criterion based on the 
SCHA in the full quantum-thermal case. As a general 
rule we can expect that the classical regime is relevant in 
the empirical range k^T/TiujQ ^1/4, which corresponds 
to TQ < 2tci, while in the opposite regime tq > 2tci 
quantum effects are dominant. 

In Fig. [1] we show the phase diagram in the full 
tq-Tci space where the instability of the SCHA occurs. 
Along the boundary line, the critical lattice fluctuations 
increase smoothly from y = 1 in the tq = case to 
?/ = 2 in the Tc\ = case. Also interesting is the depen- 
dence of the lattice fluctuations as function of Tci, namely 
the temperature (Fig. (TJ bottom panel). In the classical 
case, TQ = 0, the quadratic fluctuations y oc (u^) increase 
linearly with Td until anharmonic effects take place. An- 
harmonicity is reflected in a upturn of the temperature 
dependence of y{Tc\) and eventually in the breakdown of 
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FIG. 1: (top panel) Phase boundary of the SCHA in the rq- 
Tci space; (bottom panel) Lattice fluctuations y = A(u'^) as 
function of the classical parameter t^i for (from the bottom 
to the top) rqe = 0.2, 0.4, 0.6, . . . , 1.6, 1.8, 1.9 (we remind that 
T^Q = 2/e ). 



the solid phase for Td = 1/e and y — 1- Increasing tq 
leads not only to the presence of zero point motion quan- 
tum fluctuations at Td = 0, but also to an overall change 
of the temperature dependence of y. In particular, the 
range of the linear temperature dependence, character- 
istic of classical harmonic solids, is rapidly reduced and 
for strongly quantum solids it disappears. Lattice fluc- 
tuations are large already at T = and they are almost 
constants in a wide temperature range (note that in this 
regime anharmonic effects are in any case present due to 
quantum fluctuations) until an abrupt upturn with the 
temperature leads to the breakdown of the solid phase. 
This trends is in good qualitative agreement with re- 
cent experimental measurements^^ and Quantum Monte 
Carlo calculations We shall discuss them in details in 
Sect. IVll 
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IV. SURFACE MELTING OF QUANTUM 
SOLIDS 



After having investigated the bulk properties of quan- 
tum sohds, we analyze now the occurrence the role of 
quantum fluctuations on the surface melting. 

We can write a recursive set of equations by considering 
the quantum/thermal SCHA solution of the n-th atom 



l + 2n 



(13) 



where a;„ = y fc„ „_i „_)_i/2Af and where the local elastic 

constant fc„ „_i „_|_i is still given by Eq. ([7]). Employing 
the usual dimensionless variables tq, Tci, y™, we can thus 



write: 
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which is valid for any n > 2, while the outer atom n 
obeys the relation 



2/1 
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In order to obtain a numerical solution of Eqs. (fl^ 
(fTSj) for given tq, Tci in the stable sohd phase, we start 
by choosing a trial value of y\. The full set of {j/„} is 
thus obtained by Eqs. (|14 p - p5|) . The initial trial value 
of j/i is thus varied until yn=oo converges to its bulk value. 
Typically, this is the only physical solution, since y,i=oo 
diverges for larger values of y\ while it becomes rapidly 
negative for smaller values of y\. For tq, t^i larger than 
some critical value, the procedure does not converge for 
any value of yi, signalizing that the solid phase of the 
surface atom, described by the SCHA, is unstable. 

The resulting phase diagram, in the full tq-Tci space, 
is shown in Fig. [5] (top panel), where we compare the 
boundary of the surface melting instability (dashed line) 
with the one of the bulk melting (solid line) . For the pure 
quantum case, = 0, at zero temperature the surface in- 
stability occurs for = 0.664 where the lattice fluctu- 
ations of the outer atoms become as large as — 3.21. 
It is interesting to notice that, for 0.664 < tq < 0.736, 
the surface is unstable even at zero temperature whereas 
the bulk solid phase is always stable up to a finite tem- 
perature range. The ratio T^/Tc between the surface 
melting temperature and the temperature of bulk melt- 
ing is shown in the bottom panel of FIG. [5] showing that 
the critical temperature of surface melting can be signif- 
icantly lower than the bulk one in quantum solids. 

Before concluding this section, we would like to briefly 
compare the melting occurring at a free surface with 
other cases such as grain boundaries. In the case of a 
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bulk 
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FIG. 2: (top panel) Phase boundary for the surface melting 
instability (dashed line) compared with the bulk instability 
(solid line) in the tq-Tci space; (bottom panel) Ratio between 
surface melting temperature T^"^ and bulk melting temper- 
ature Tc as function of tq. For rge even the bulk phase is 
unstable. For 1.12 < rqe < 2 the system is in a quantum 
regime where UtuT^^ /fiuJo < 1/4. 



free surface, in going from Eq. p4)l to Eq. (|T5|) , we have 
dropped in Eq. (jlSp the contribution of the n — Q atom. 
We note that the same results would be obtained in Eq. 

considering n = I and assuming the lattice fluctua- 
tions at the site n = to be infinite, namely yn=o = oo. 
This latter condition would be obtained by the harmonic 
oscillator solution of at the site n = with a van- 
ishing clastic constant kn^n-i,n+i, and it express nothing 
more than the condition that atoms for n < \ are not in 
a solid arrangement but in a gaseous phase. 

An intermediate situation is encountered when melt- 
ing at grain boundary interfaces is considered. In this 
case the outer atom n = 1 of a grain would not probe a 
free surface at the site n ~ Q, but it will interact with a 
lattice environment with a different arrangement. These 
two situations can be described by a similar set of re- 
cursion relations (|14p but with different boundary con- 
ditions: in the free surface case boundary conditions at 
site n = will be described by a completely soft oscillator 
kn,n-i,n+i = 0, signalizing that bulk solid is interfaced 
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with a free gaseous phase; on the other hand, in the case 
of grain boundaries, the outer atom n = 1 wiU stiU probe 
a crystal structure for n < 0, although with a different 
arrangement. The boundary conditions at site rt = 
will be still described thus by Eq. (fT3)l . but with a not 
completely soft mode. We expect thus that melting pro- 
cesses occur as well at grain boundaries as in the case of 
free surface. From the mathematical point of view, this 
situation is identical to the case of quantum isotopic sub- 
stitutions, and it will be discussed in details in the next 
section. 



QUANTUM MELTING DRIVEN BY 
ISOTOPIC IMPURITIES 



In this section we address the problem of the solid 
phase stability close to a single local isotopic substitution 
embedded in a perfect lattice structure. In the SCHA ap- 
proach, local stability of the solid phase is given by the 
solution of Eq. It is easy to check that, in the clas- 

sical limit k^T ^ fii^n, the dependence on the atomic 
mass M in Eq. (fTH]) drops out, so that different isotope 
solids should probe the same stability conditions. On the 
other hand, the mere observation of a different melting 
line for ''He and •^He is a direct evidence that helium is 
in a quantum regimei^ii^ Different isotopes are thus ex- 
pected to affect the bulk solid phase stability. We expect 
the same at the local level. 

In the following we shall consider the case of a isolate 
substitution with a lighter isotope in a host matrix of 
heavier atoms. Quantum fluctuations in the two cases 
will be ruled locally by the parameters tl — Xfi/2y/koMh 
th = Xh/2^/koMB_, respectively for the lighter (L) and for 
the heavier (H) atoms. To study the stability of the solid 
phase close to this isotopic quantum impurity, we can still 
employ the recursive relations pij) . namely for n < —1, 
n > 1 we set tq = th, whereas for n = (quantum 
isotope impurity) we have tq = tl- We shall consider 
the representative case of a ^ He imp urity embedded in 
''He solid. In this case tl/tr = a/4/3. 

In Fig. [ni^top panel) we show the phase diagram of 
the lattice instability of the host ^He solid close to the 
quantum isotopic "^He impurity. It is instructive to com- 
pare the classical limit tq = with the pure quantum 
one Tci = 0. In the first case lattice fluctuations of the 
guest atom, as well as of the host atoms, are indepen- 
dent on the relative atomic mass and they depend only 
on the temperature. As a consequence, the solid phase 
close to the guest atom is completely unaffected by the 
isotopic substitution. A quite different situation occurs 
in the highly quantum regime Td = 0. In this case lo- 
cal quantum lattice fluctuations of the lighter guest atom 
can be significantly enhanced due to its reduced atomic 
mass, and they can be sufficiently large to induce a local 
melting of the host solid phase. At Td — this occurs 
for Th > 0.681, not much higher than in the case of a 
free surface truncation (tq > 0.664). Note that Fig. [3] 
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FIG. 3: (top panel) Pliase boundary for tlie lattice instability 
around a quantum isotopic substitution witli tl/th = \/ 4/3 
(dashed line) compared with the bulk instability (solid line); 
(bottom panel) Ratio between melting lattice temperature 
r^'™P around the quantum impurity and bulk melting tem- 
perature Tc as function of host quantum parameter th- In the 
quantum regime 1.39 < rne < 2, where ksT^'^'^'^ /hujo < 1/4, 
the local melting temperature around the quantum impurity 
is sensible lower than in the bulk, and for 1.85 < rne < 2 solid 
phase around isotopic quantum impurities is unstable even at 
T = 0. 



defines a region (quantum impurity melting) where solid 
phase is still stable in the bulk but local quantum lat- 
tice fluctuations break down the solid phase close to the 
isotopic substitution. On the physical ground we can ex- 
pect liquid bubbles of host atoms to appear close to the 
guest isotope. Unfortunately, since the present analy- 
sis is only related to the stability condition of the solid 
phase, we are not able to estimate the size of the liquid 
bubble, and more sophisticated approaches are needed. 
It is interesting to note that, for quantum solids, the crit- 
ical temperature 7)9"°? for the local stability of the solid 
phase close to the quantum isotope impurity is reduced 
with respect to the bulk Tc- This is shown in the bottom 
panel of Fig. [3] where the ratio between the local Tc close 
to the impurity and the bulk Tc is plotted as function of 
the quantum degree of the system, parametrized by th- 
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In the quantum regime, where T^™? < the lo- 

cal melting temperature T^™p can be significantly lower 
than the one in the bulk Tc, and, for 1.39 < thc < 2, 
we expect a quantum isotopic impurity to induce local 
melting down to T = 0, although the bulk phase is still 
stable. 



VI. DISCUSSION AND CONCLUSIONS 

In this paper we have investigated the stability of quan- 
tum solids with respect to surface melting and to isotopic 
quantum substitutions. Both these phenomena can be es- 
sentially related to the amount of lattice fluctuations, and 
they can be driven thus by thermal fluctuations as well as 
by the zero point quantum motion. We have shown that 
the effects of isotopic impurities and surface melting are 
strongly enhanced in quantum solids. In particular we 
show that when quantum fluctuations are dominant in 
quantum solids the solid phase can be rapidly destroyed 
on the surface and close to quantum impurities at tem- 
peratures much smaller than for the bulk melting. 

Helium solids are the natural candidates where the 
quantum instabilities of surface or interface can occur. 
The actual relevance of these quantum melting effects 
are of course ruled by the magnitude of the quantum lat- 
tice fluctuations which are parametrized in our model by 
the quantity tq. An accurate calculation of the quan- 
tum lattice fluctuations as a function of the tempera- 
ture in '*He and ^He solids has been provided recently, 
by using of Quantum Monte Carlo (QMC) techniques, 
by Draeger and Ceperley in Ref. HO, in excellent agree- 
ment with the experimental datai^ Quite interestingly, 
they find that the mean square lattice displacement {v?)t 
does not follow at low temperature an harmonic behavior 
{v?)t — {u'^)t=o + OiT^, but rather a more shallow one 

Ref. |20| represents a suitable source to estimate an ef- 
fective value of Tq representative of solid helium. To this 
aim we fit the temperature dependence of the QMC data 
of Ref. 113 with our quantum SCHA model described by 
Eq. (|12p . where only two independent fitting parameters 
appear, namely A and ko (remind that Tc\ = Xk^T/kQ, 
TQ = Xh/2y/k^). The fit of our quantum SCHA [Eq. 
P^ ] compared with the QMC data is shown in Fig. |3]for 
three representative cases where the number of numeri- 
cal data is larger than the number of independent fitting 
parameters to guarantee the significance of the fitting 
procedure. Also shown is the fit with a purely harmonic 
model obtained by setting A = 0. The extracted values 
of A and fco , as well as of the corresponding tq and of the 
anharmonic renormalized phonon frequency at T = cDo 
are reported in Table [H where also we report the critical 
temperature Tc for the solid phase bulk instability eval- 
uated within the SCHA and the experimental melting 
temperature T^Pfiii^^ 

It is worth to comment about the temperature behav- 
ior of the QMC data compared with the harmonic (A = 0) 
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FIG. 4: Lattice fluctuations (u^) evaluated within the SCHA 
(solid lines) as function of temperature for different helium 
solid conditions compared with Quantum Monte Carlo data 
of Ref. [20I . Values of fco and A in SCHA obtained by fitting 
QMC data are reported in Table |T] Also shown are the purely 
harmonic fitting of the QMC data with a Einstein and a Debye 
model. 



and anharmonic SCHA fit. An important point to be 
here underlined is that QMC results show a large mean 
square lattice displacement at zero temperature all to- 
gether with a rapidly turn up of {v?) close to the solid 
bulk instability. As we have discussed in Sect. IIIIl this 
is a characteristic trend of highly quantum solids. On 
the other hand, this behavior is poorly reproduced by 
a purely harmonic model where the amount of the lat- 





hep ''He 


fee ■'He 


fee =^He 


Vo (cm^/mole) 


12.12 


10.98 


11.54 


fco (meV/A^) 


110 ± 10 


140 ± 10 


150 ± 10 


A (A-) 


14 ± 1 


15.2 ± 0.8 


14.7 ±0.7 


TQ 


0.69 ±0.08 


0.66 ±0.06 


0.70 ± 0.06 


Cjo (meV) 


4.6 ±0.1 


5.4 ±0.4 


6.1 ±0.5 


Tc (K) 


14 ±4 


20 ±4 


18 ±4 


T-P (K) 


~ 15 


~ 21 


~ 22 



TABLE I: Values of fcg and A in SCHA obtained by fitting the 
QMC data of Ref. 20 for three representative helium solids, 
namely: hep ^He at molar volume Vo — 12.12 cm^/mole, 
fee *He at molar volume Vb = 10.98 cm'^/mole, and fee ^He 
at molar volume Vo = 11.54 cm^/mole. Also reported are 
the corresponding values of tq, the renormalized phonon fre- 
quency Cbo and the predicted critical temperature Tc of the 
solid phase bulk instability compared with the experimental 
melting temperature 21^ 
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tice fluctuations at T = is inversely proportional to 
the temperature dependence. This is even more true if a 
Debye model would be employed since the temperature 
dependence of a Debye model is even more shallow than 
in the Einstein case. 

The strong quantum degree of solid helium, qualita- 
tively predicted by these arguments, is confirmed by the 
numerical analysis of the SCHA fit which predicts a quan- 
tum parameter tq in the range tq ~ 0.66 — 0.7 for the 
three samples here considered. The robustness of our fits 
is confirmed by the nice agreement between the critical 
temperature for the bulk instability of the solid phase 
estimated by the SCHA and the experimental melting 
temperature. 

These results have important consequences with re- 
spect to the surface/grain-boundary melting instability 
and local melting induced by quantum isotopic impuri- 
ties. The values of tq ~ 0.69, for the low pressure/ high 
molar volume Vq — 12.12 cm'^/mole, is safely larger than 
the value Tq^ ~ 0.664 where surface melting occurs at 
zero temperature, and also or the same order and slightly 
larger even than Tq^ ~ 0.681 where isotopic impurity 
induced melting also occurs at zero temperature. Al- 
though these estimates have to be meant only indicative 
of the quantum degree of helium solid, they clearly point 
out that quantum anharmonic effects are large enough 
in solid helium, for these or larger molar volumes, to en- 
force surface melting and local melting close to quantum 
impurities down to zero temperature. Quantum Monte 
Carlo simulations have actually confirmed premelting at 



surface between helium solid and Vycor walls-^ and in- 
ternal interfaces of a pure helium system^^, although not 
all possible interfaces undergoes a solid/liquid transition. 

These results shed an interesting light also on the 
recent report of the Non-Classical Rotational Iner- 
tia (NCRI) observed in ^He-^iSl While it was ini- 
tially claimed to be an evidence of a supersolid (SS) 
phase, subsequent experiments showed a strong de- 
pendence of the NCRI on the annealing process^^ on 
the presence of grain boundaries;^ on the amount of 
■^He concentration2Zi^i2i as well as on the freezing 
procedurei^i^ These observations give rise to an alter- 
native hypothesis to the SS phase, namely, that a liquid 
phase is confined at the grain boundaries and that mass 
flow is related to superfluidity of the liquid component 
Our results confirm this scenario and shed new perspec- 
tives about the role of disorder/grain boundaries in solid 
helium. In particular we provide a natural explanation 
for the existence of a liquid (and thus probably super- 
fluid) phase at the grain boundaries and we predict a local 
liquid phase also around "^He impurities. Local melting 
close to isotopic '^He impurities should be thus explicitly 
considered. 
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